Second Thresholds in BEC-BCS-Laser Crossover of Exciton-Polariton Systems 
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The mechanism of second thresholds observed in several experiments is theoretically revealed by 
studying the BEC-BCS-laser crossover in exciton-polariton systems. We found that there are two 
different types for the second thresholds; one is a crossover within quasi-equilibrium phases and the 
other is into non-equilibrium (lasing). In both cases, the light-induced band renormalization causes 
gaps in the conduction and valence bands, which indicates the exsistence of bound electron-hole 
pairs in contrast to earlier expectations. We also show that these two types can be distinguished by 
the gain spectra. 
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In semiconductor exciton-polariton systems, Bosc- 
Einstein condensation (BEC) of exciton-polaritons has 
been observed in recent years [D-Q- A hot issue is, now, 
how the exciton-polariton BEC, a thermal equilibrium 
phenomenon, changes into the lasing operation resulting 
from the electron-hole (e-h) plasma gain [fl, which is es- 
sentially a non-equilibrium phenomenon 0, 3 . Earlier 
experiments show that there are two distinct thresholds 
when increasing the excitation density: the first thresh- 
old is the critical density for the BEC {jj [l(| and the sec- 
ond one is recognized as the standard lasing 
most cases, the second-threshold mechanism is explained 
by a shift into the weak coupling regime due to dissoci- 
ations of bound e-h pairs into the e-h plasma. However, 
there is no convincing discussion why such dissociations 
lead to non-equilibriation of the system essential for las- 
ing. As another possibility, a new ordered state involving 
Bardeen- Cooper- Schrieffer (BCS) -like correlation is also 
speculated [15(. The second threshold is, thus, currently 
subject to intense debate. In this letter, our purpose is to 
reveal the mechanism of the second threshold by study- 
ing the BEC-BCS-laser crossover theories [l6-18|. As a 
result, we found that there are two different types for the 
second threshold; one is a crossover into photonic polari- 
ton BEC (quasi-equilibrium phase) [3, [2(J and the other 
is into lasing (non-equilibrium). In both cases, the light- 
induced band renormalization causes gaps inside the con- 
duction and valence bands, which indicates that there are 
still light-induced e-h pairs even after the second thresh- 
olds, in contrast to the above scenario. We also show that 
these two types can be distinguished by the gain spectra. 

Our starting Hamiltonian is H = Hg + + H^r, 
where 
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are the system, reservoir, and their interaction Hamil- 
tonians, respectively (HI- Here, a, a' € {c, v}, H.c. 
is the hermitian conjugate, c c .fc (Cv,fe) is the conduc- 
tion (valence) band electron annihilation operators with 
wavenumber fe, U' q is the Coulomb interaction, g is the 
light-matter coupling constant, £ c / v .fc = ^> c /v,k T H fi/2 
is the electronic dispersion measured from ±h~ 1 n/2, a q 
is the cavity photon annihilation operator with wavenum- 
ber q, £ p h,q = ^pii.q — fr ft is the photonic disper- 
sion measured from h fi, and H fj, is an oscillation 
frequency of the photon and polarization fields, which 
corresponds to the energy of the main peak in photolu- 
minescence. Similarly, b c k and b v ^ k denote fermion an- 
nihilation operators of pumping baths and \P P is a bo- 
son annihilation operator of free-space vacuum fields. In 
this model, Eqs. ([2]) and ((3]) are responsible for the in- 
coherent fermionic pumping and photon decay. Based 
on the above Hamiltonians, we focus on steady states 
described by the polarization function p k = (c^ fc c Cj fc), 

the number of electrons n c ^ k = (cj k c c _ k ) and holes 

nh,-fc = 1 — (cj fcC v .fc) , and a coherent photon field formed 
in the q = state (a q ) = Sq^ao with the oscillation fre- 
quency hT 1 [i. These are the minimum variables for de- 
scribing the BEC, BCS states, and semiconductor lasers. 

The band renormalization of the e-h system can, 
then, be conveniently studied by the poles of the single- 
particle spectral function A aa i [y; k) = i(G^ a , iv; k) — 
G< a ,(v-k)), where iHG> a ,(tt';k) = (c Q , fc (t)^, fc (t')> and 

-iSG^,(tf';fc) = (c^, )fe (t')ca,fe(i)) are the greater and 
lesser Green's functions in the time domain, respectively. 
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Within the Hartree-Fock (HF) approximation, a stan- 
dard Green's function technique yields 



/, B (-v) 



A cv {v;k) = 2u* k v k [L(v,E k )-L(v,-E k ) 



(4) 



A cc/vv {v; k) = 2\u k \ 2 L(v, T £ fe ) + 2\v k \ 2 L{v, ±E k ), (5) 

where y/2u k = [1 + £+ k /E k }^ 2 and V2v k = e^[l - 

&,J E -] 1/a with E * = [(ei lh ) a + l^| 2 ] 1/2 are the Bo- 
goliubov coefficients, L(v,±E k ) = h^/[{hv — h4~ ilk ± 
hE k ) 2 + (Try) 2 ] is the Lorentz function, and 7 is the ther- 
malization rate of the e-h system [2l| . Here, £ e /h,fe = 
^c/h.fc + - ft~V2 with |± fe = (| e ,fe ± | h)fc )/2 

describes the single particle energy renormalized by the 



Coulomb interactions hE^J^, = 



k n c/h,k> 



and 



A k is a composite order parameter defined as A k = 
\A k \e i9k = g*a„+h- 1 £ fe , U k ,_ k p k >. We note that, in the 
derivation, the notation is transformed into the e-h pic- 
ture by uj Ctk = u Ctk and w h ,fc = -u y ,k + ^ J2k' U L'-k- 
In Eqs. ([3]) and ([5]), there are remarkable similarities 
to superconductivities [22j. This is because the HF ap- 
proximation is essentially the same as the mean-field ap- 
proximation for the system Hamiltonian (Eq. (HJ), which 
can be diagonalizcd by the Bogoliubov transformations. 
It is then clear that min[2?i£'fe] represents the gap energy 
opened at fj,/2 in the renormalized conduction and va- 
lence bands (typically Fig.Q](a) and (b)). Such a picture 
is well-known for e.g. the BEC-BCS crossover physics 
but, now, one should notice that Eqs. (|4]) and ([5]) are 
also applicable for lasing because thermal equilibrium is 
not required for their derivations. In the case of lasing, (i 
is not the chemical potential but the laser frequency and 
the gap is opened whenever lasing because ao ^ and 
p k result in mm[2HE k ] 7^ 0. The origin of the gap is 
analogous to the Rabi splitting in resonance fluorescence 
[23M25I ] . This can be understood from the expression of 
mm[2HEk] = 2h\g\y/nph {n p h = |ao| 2 ) obtained by as- 
suming free electrons (U' q = 0) with fx > E g (E g ; the 
bare band gap energy), which is equivalent to the Rabi 
frequency in resonance fluorescence. Furethermore, the 
existence of the gap indicates that light-induced e-h pairs 
do exist whenever lasing even though there is no e-h pair 
before lasing. This is one of our important results despite 
the quite simple analysis. For later convenience, two typ- 
ical situations for large and small gap energies are shown 
in Fig.[T](a) and (b), respectively. In Fig.[T](a), the renor- 
malized conduction band has a gap around k ~ with 
flattened dispersions because e-h band mixing occurs for 
large fc-regions. In contrast, in Fig. [T](b), the renormal- 
ization is mainly focused on particular fc-regions. In both 
cases, the renormalized bands have gaps at \xj2 and the 
same holds for the valence band (not shown). 

In order to discuss the second-threshold mechanism, 
however, the unknown variables A k (ao and p k ), n Cik , 
rth,fc, and \i in Eqs. (@| and (JSJ) should be determined 
in a comprehensive way including BEC, BCS and laser 




FIG. 1: Renormalized conduction bands for (a) large (\fj, — 
E*\ < min[2ft_E fe ]) and (b) small (mm[2hE k ] < fj, — E*) gap 
energies. Here, E* = E g + hEf k = + hE^^o an d the gray 
dashed lines show the energies by ignoring A k in E k . Rela- 
tions to /^ h (±i/) are also illustrated in panel (c) for quasi- 
equilibrium (min[2hE k ] > /ib — n + 2h"/-\-2kBT) and panel (d) 
for lasing (/j,b — fJ. > min[2hE k ] + 2fry + 2k,B.T). ui a ^ = Wh,k 
and /ij? = fj,® are assumed. 

physics flil- 18|. Within the HF approximation, the si- 
multaneous equations can be written as 



hA k 
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where U k ,' k = /%| 2 /(£ P h,o - i«) + U' k ,_ k is the effective 
e-h attraction and k is the photon loss rate [2l| . Eaxt ex- 
pressions of /^hffe S (^) and f^ s k (v) are given in supple- 
mentary [2l[ , which are described by the bath Fermi dis- 
tribution functions /^ h (z^) = [exp(P(hv—(n^ h — fj,/2))) + 
I] -1 with the inverse temperature /3 = and the 

chemical potentials £t B / h - The important point here is 
that, by assuming u c ^ k = Wh,fe and a charge neutrality 

l"? = Mh with Mb = + Mh> E q s - © and © can 
recover the BCS gap equation when (I) min[2hE k ] > 
Mb~ /-i+2h r f+2kBT (quasi-equilibrium) , while there are k- 
regions described by the Maxwell-Semiconductor-Bloch 
equations when (II) /ib — M ^ min[2hE k ] + 2thy + 2k^T 
(lasing states; BCS-coupled lasing in a strict sense). The 
physical mea ning s of these conditions are discussed in 
detail in Ref. [I8[ and not repeated here. Instead, these 
conditions are illustrated in Fig. Q] (c) and (d) in relation 
to /3 ll (±i / ) for quasi-equilibrium and lasing conditions, 
respectively. It is, then, clear that the system enters into 
lasing phases when — M roughly goes beyond the en- 
ergy gap of mm^2hE k ] by ignoring the broadening due to 
7 and T. 

Based on the above formalism, we have performed nu- 
merical calculations, where the cavity level (= ?iw p h,o) 



3 




Binding energy 

of polariton 

10-' 10° 10 1 10 2 
Pumping :/u B - E LP (meV) 

FIG. 2: Numerical solutions of (a) the number of photons 
|cto | 2 , (b) the oscillation frequency /x, and (c) the gap energy 
min[2/LEfc] as a function /xb for Tik = 0.1 ixeV and 100 /ieV. 
Red and blue colors are used when satisfying the condition (I) 
for quasi-equilibrium and the condition (II) for lasing in the 
text, respectively. Green colors are used when both of them 
are not satisfied. For comparison, black solid lines show the 
results by thermal-equilibrium theories 0, [2(| . 

is in resonance with the (IS) exciton level located at 
10 meV below E g and the lower polariton level E'lp is 
formed at 20 meV below E g |2lj . For Hk, we have used 
values of 0.1 /xeV and 100 /icV to study the effects of non- 
equilibrium. We note, however, that Hk = 100 /xeV is a 
reasonable value in current experiments. Fig.[5]shows the 
calculated results of |ao| 2 , p, and mm[2HEk] as a function 
of fiB, the pumping parameter. In the case of the equilib- 
rium theories, |ao| 2 diverges in the limit of /xb — ► feph.o 
because it is preferable to increase photons rather than 
electrons and holes due to the phase space filling effects. 
As a result, the photonic polariton BEC is achieved by 
the photon-mediated e-h attraction 0, [20|. In contrast, 
in the case with finite pumping and losses (plots), the 
behaviors are different in many aspects. Focusing on the 
plots for Hk = 0.1 /xeV, two distinct thresholds can be 
seen (p B - E LP w 5.0 x 10" 1 meV and 6.0 x 10° meV) 
in Fig. 03 (a). At the same time, fx is gradually blue- 
shifted from _Elp and then approaches the bare cavity 
resonance. Similar qualitative behavior also can be seen 
for Hk = 100 /xeV. These behaviors are consistent with 
experiments. However, there is a crucial difference be- 



1.0 

e 

■B 0.5 

C3 
N 

'§ o.o 

"o 1.0 

a, 

-a 

S 0.5 



e 

■c 



0.0 
1.0 



0.5 



0.0 



hK= 0.1 u eV 


hK= 100 ueV 


: (a) 

■ >i 


: (d) 


' (b) 




(c) 

" .... 1 ... . 





► Before 
second threshold 



After 

second threshold 



0.0 



0.5 



0.0 



0.5 



1.0 



Wavenumber k ( X 10 9 nr 1 ) 

FIG. 3: Calculated distributions n c ,fe and polarizations \pk\ 
corresponding to the triangles indicated by a-f in Fig. [2(a). 
Panel (a)-(c) are for Hk — 0.1 /xeV and panel (d)-(f) are for 
Tin = 100 /xeV. Arrows show the kinetic hole burning. 

tween the two; according to the above-mentioned condi- 
tions (I) and (II), all plots are in quasi-equilibrium for 
Hk = 0.1 fieV but there are plots (blue) in lasing for 
Hk = 100 fxeY after the second threshold. 

The difference is also reflected in n c ^ and pu, as shown 
in Fig. [3l Before the second thresholds, n Cj fe and Pk 
for Hk = 0.1 /LteV (Fig. [3] (a)) are similar to those for 
Hk = 100 /j,eV (Fig. [3] (d)). However, after the second 
thresholds, n c ^k and pk are quite different, depending on 
the value of Hk. In the case of Hk — 0.1 /icV, n c .fe mono- 
tonically decreases as a function of k and pk has a plateau 
w 0.5 (Fig. [3] (b) and (c)), which are the same features 
as the photonic polariton BEC in thermal equilibrium 
[3 In contrast, in the case of Hk = 100 /xeV, the 
kinetic hole burning appears as a signature of lasing and 
the Fermi surface is formed with the population inversion 
?i c ,fe > 0.5 (Fig. 03(e)). For larger /xb, such behaviors be- 
come much more pronounced (Fig. 03(f))- These results 
directly show that the second thresholds for Hk — 0.1 /xeV 
and 100 /xeV in Fig. 03 (a) are formed by different mech- 
anisms. 

In fact, for Hk = 0.1 peV, the second threshold is 
formed by the same mechanism as the photon diver- 
gence in the equilibrium theories, and therefore, it results 
from the crossover into the photonic polariton BEC. In 
the present case, there are finite losses of cavity photons 
even if the system is in quasi-equilibrium. As a result, 
the divergence is avoided and the second threshold ap- 
pears instead. After the second threshold, the monotonic 
increase of min[2ft£'fc] (Fig. 03(c)) indicates the enhance- 
ment of the light-induced e-h paring and expands the flat- 
tened region of dispersion in Fig.Q](a). It is, then, clear 
that the plateau of pk ~ 0.5 in Fig. 03(b) and (c) is formed 
by the e-h mixing around such flattened dispersions. In 
the case of Hk = 100 /xeV, on the other hand, the second 
threshold is related to the crossover into lasing, explained 
as follows. Before the second threshold, the system stays 
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in quasi-equilibrium (red circles in Fig. [2]) , where the re- 
lationship between the rcnormalizcd band and the pump- 
ing baths is well expressed in Fig.Q](c). In this situation, 
the pumping is blocked inside the gap min[2hE k ]. How- 
ever, by increasing the pumping /jb, Mb — A* exceed the 
gap, — H min[2?LEfc], and then, electrons above the 
gap can be supplied suddenly. Such a feeding mecha- 
nism causes a rapid increase of photons, resulting in the 
second threshold. Here, by ignoring the effects of 7 and 
T, this situation Mb ~ M ^ mm[2hEk] is equivalent to 
the above-described condition (II) for the lasing phases. 
Consequently, the second threshold is accompanied by 
the change into lasing (non-equilibrium). By increasing 
the pumping further, fi is fixed around the cavity (Fig. [2] 
(b)), min[2HEk] is decreased (Fig. [2] (c)), and the effec- 
tive band gap E* = E g + hE*g* + hS^% shrinks, 
of course. The lasing situation is then well captured in 
Fig.[T](d), where the gap min[2/lEfc] is decreased but still 
opened around the laser frequency. The decrease of the 
gap for |Ub — A* ~ mm[2hEk] implies that the particle 
flux hb — M beyond min[2/lEfc] act toward e-h pair break- 
ing but the e-h pairs cannot be fully dissociated because 
mm[2HEk] 7^ 0. As a result, light-induced e-h pairs are 
still formed around the laser frequency, typically around 
the energy regions of the kinetic hole burning (Fig. [3] (c) 
and (f)). This is, in turn, somewhat analogous to the 
e-h Cooper pairs formed around the Fermi energy, i.e. 
weakly correlated e-h pairs in momentum space. The 
difference is that the e-h pairs are formed around the 
laser frequency rather than the Fermi energy. 

These results indicate that it would be reasonable to 
explain the second thresholds reported in current exper- 
iments by the crossover into lasing because k = 100 A*eV 
is a reasonable value for them, in agreement with ear- 
lier explanations [H, [H-14|. Our theory, however, show 
that the crossover is not accompanied by the dissocia- 
tions of bound e-h pairs. Instead, the pairing mechanism 
changes into the light-induced one around the laser fre- 
quency. This is in contrast to the commonly accepted 
ideas but a natural picture of lasing. 

We have thus discussed the two different types of the 
second threshold. However, it is difficult to directly dis- 
tinguish them by the excitation dependence of the num- 
ber of photons (Fig. [2] (a)), in principle. Therefore, we 
finally study the measurable optical gain spectra G(iS) 
26L 27j by assuming an additional perturbative Hamilto- 



te=0.1|jeV 



hK= 100 neV 



nian H'(t) = -F(t) £^(2^,* + H.C.). Here, F(t) 
is the weak light field irradiated from the outside and 
d cv is the dipole matrix element. Within the linear re- 
sponse |28j], G{uj) is estimated with the ladder approxi- 
mation [29(. Figure H] (a)-(f) shows the gain spectra cor- 
responding to Fig. [3] (a)-(f), respectively. In the case of 
Hk = 0.1 /^eV, two absorption peaks can be found, which 
result from the two flattened dispersions shown in Fig. [1] 
(a). Therefore, the separation of the peaks corresponds 
to min[4S£'fc], the sum of the gaps in the conduction and 
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FIG. 4: Optical gain spectra. Panel (a)-(f) corresponds to 
Fig- E] (a)-(f), respectively. Positive (red) and negative (blue) 
values represents the gain and absorption, respectively. 

valence bands. Here, we note that absorption dominates 
the spectra because there is no or little population inver- 
sion (n e ,fc > 0.5) for the condensed phases in equilibrium 
(Fig. [3] (a)-(c)). In the case of Hk = 100 fieV, however, 
gain appears when the system enters into the lasing phase 
(Fig. 2] (e) and (f)) although absorption still dominates 
in the quasi-equilibrium case (Fig. U (d)). The spectral 
hole (or gap) with a separation of mm{AKEk\ in Fig.|4](f) 
reflects the gap formed in the renormalized band (Fig. Q] 
(d)). The existence of the gain after the second threshold 
is due to the population inversion for lasing (Fig. |3] (c) 
and (f)) and, as a result, can be used to distinguish the 
two types of the second threshold in experiments. 

To summarize, we have shown that there are two dif- 
ferent types for the second threshold. In both cases, dis- 
sociations of bound e-h pairs do not occur due to the 
light-induced pairing, in contrast to earlier expectations. 
The gain spectra are also studied and the existence of 
the gain would be useful to distinguish the two different 
types of the second threshold. 
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Supplemental material 

In this supplementary material, we describe a few details of our formalism in order to make the paper self-contained. 
We also describe the parameters in our numerical calculations and compare the results with our own previous study 
[Yamaguchi et al, New J. Phys. 14, 065001 (2012)]. 



A few details of the formalism 



In the main text, the thermalization rate and the photon loss rate are described by 7 and k, which result from the 
system-bath interaction Hamiltonian Eq. ([3]). With the definitions of the density of states 

D*(u>)=J2 S «*-^> J&M = J>( W P- W ). (8) 

k p 

7 and k have the following relations with and Cq : 

7 = la,k = nlr^D^u), K = Kq = ^|C q | 2 ^ h (^). (9) 

Here, the dependence on the wavenumber is neglected for simplicity. 

In addition, in Eqs. (j6)) and (J7J, the distributions of /^j? S (^) an d /<!/h^f M are defined as 

flTM = k/» - 4 B (-^)] + + / B (-*) - 1]^^, (io) 

2 2 v- £ ehfc 

/Xf M = /c/hW'fc/h^M + [1 - / h B / e (-^)][l - VMM], (11) 

where r/e/h.feC 1 ') is 



S (, + 6Ve, fc ) 2 + 7 2 + W' (12) 



(^ + £h/e,fc) 2 +7 2 

h |h/e,fe) 2 + T 2 + I- 

with < f] c /h.k( l/ ) < 1- Then, from Eqs. (jSJ), (JTJ), ([TU|) and (|11[) . it can be found that the band renormalization 
effects arc included not only in A aa i(v;k) but also in f^nf'^iv). For example, in Eq. (|TT|) . T) e /b,k( 1 ') = 1 an d 
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/c/h S fc S (^) = flifiSy) when Z\fc = 0. This means that electrons and holes are thermalized into the simple Fermi 
distribution. However, there exists e-h band mixing when 7^ 0, and therefore f^/huty) * s influenced by 1 — f^/J^v) 
with a weighting factor of 1 — rjg/^ ^u), as seen in Eq. (JXTJ) . 



Parameters in the numerical calculations 



In our numerical calculations, the fc-dependence of Z\fc is eliminated by using a contact potential U' q = U and the 
other parameters are fiu; ,fc = fiwh.fc = h 2 k 2 /2m + E g /2, m — 0.068mo (mo is the free electron mass), zZg = /i^ , T = 
10 K, and fvy = 4 meV. In this context, we note that the calculations are qualitative even though these parameters 
are taken as realistic as possible. Here, the contact potential and the coupling constant are, respectively, set as 
U = 2.66 x 1CT 10 eV and hg = 6.29 x 1CT 7 eV with cut-off wavenumber fc c = 1.36 x 10 9 rrT 1 . In this case, the (IS) 
exciton level (= E ex ) is formed at 10 meV below E g (E CK = E g — 10 meV) and the lower polariton level E^p is created 
at 20 meV below E g (E^p = E s — 20 meV) under the resonant condition fiw p h,o = E ex , as described in the main text. 



Comparisons with our previous results 



In the present study, there are mainly two differences from the previous conditions of numerical calculations reported 
in Ref. [Sl|. The first difference is the cavity detuning; the cavity resonance foj p h,o is tuned to the exciton level E^ 
in the present case (fiw p h,o = E cx = E g — 10 meV), whereas it is largely detuned from the exciton level in Ref. |Sl( 
(fiw p h,o = E cx + 40 meV = E g + 30 meV). The second difference is the value of the coupling constant hg which 
determines the position of the lower polariton level under the resonant condition; the lower polariton level is located 
at 10 meV below the exction level in the present case (-Elp = E cx — 10 meV), while it is at 5 meV below the exciton 
level in Ref. fs~0 ] (£bp = E CK — 5 meV). Based on these differences S2|. we compare the previous and present numerical 
results here. 

In the previous case, as a result of the large detuning and the small coupling constant, the cavity has little influence 
on the behavior of excitons formed in the low density regime. Only in a high density regime, the cavity has large 
influence on the behavior of the condensation but excitons no longer exist in the high density regime. Such signatures 



can be found in Fig. 10 and Fig. 11 of Ref. [Si] . Red and blue plots in Fig. 10 (a) show that the order parameter 
HA comes up around /ib — E g = 10 meV when /zb is increased, as indicated by the arrow (A). This means that the 
condensation starts when /zb reaches around the exciton level, i.e. p,B = E cx = E g — 10 meV and there is, indeed, 
little influence of the cavity. In addition, Fig. 11 (A) shows that n CJ t and \p^\ monotonically decrease in the fc-space. 
Therefore, the condensation indicated by (A) in Fig. 10 (a) is identified as the exciton BEC. The negligibly small 
photonic fraction around the point (A) in Fig. 11 is also consistent with this interpretation. Here, we note that the 
photonic fraction is still small even though /iB is increased up to around the point (B). Hence, the effect of the cavity 
is still small at this point. However, in contrast to Fig. 11 (A), Fig. 11 (B) shows that n ^ has a rounded shape 
of the Fermi function and \pk | has a peak around the Fermi surface. Therefore, the condensation around the point 
(B) can be classified as the electron-hole (e-h) BCS phase. It should be noted that, at this stage, there is no exciton 
anymore because the density is increased as much as the Fermi surface is formed. This means that the cavity has 
little influence on the crossover from the exciton BEC to the e-h BCS phase. In this regime, it is natural that there 
is little effect of the cavity loss k in Fig. 10 (a)-(c) (the difference between the red and blue plots for /zb < ftw P h,o)- 
This effect becomes apparent only after /xb is increased up to around the cavity resonance (/xb > Sw P h,o) but excitons 
no longer exist in this situation, as described above. 

In the present case, on the other hand, the cavity has large influence on the condensation even in the low density 
regime due to the resonance condition and the larger coupling constant. In Fig. [2] (a), it can be found that the 
condensation occurs when /ib is around £"lp rather than the exciton level. Given the monotonic decrease of n c _k 
and \pk\ in fc-space, as shown in Fig. |3](a) and (d), the condensation can be identified as the exciton-polariton BEC. 
As a result, the cavity loss k also has a great impact on the behavior of the condensation, as evidenced in Fig. [5] 
(a)-(c) (the differences between the plots of n = 0.1 /xeV and k = 100 //eV). This is in contrast to the previous results 
where k plays no role in the low density regime (pb ^ ^ P h,o in Fig. 10 (a)-(c) of Ref. [sTJ] ) . However, when /j,b 
is largely increased, the situation is not so different from the previous calculations. This is because the band gap 
renormalization decreases the effective band gap E* and the cavity level enters into the conduction and valence bands. 
Therefore, for example, the behaviors of n e ,fc and \pk\ after the second thresholds in Fig. [3] arc similar to those in 
Fig. 11 (E) and (H) of Ref. |Slj . Thus, from the comparison of the present and previous numerical results, we can 
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learn that the effect of the cavity in the low density regime becomes more important in the present case than in the 
previous one. 



* E-mail: yamaguchi@acty.phys.sci.osaka-u.ac.jp 
[SI] M. Yamaguchi, K. Kamide, T. Ogawa, and Y. Yamamoto, New J. Phys. 14, 065001 (2012). 

[S2] T = 10 K is also different from the previous condition T = K but this difference does not play any significant role because 
&bT is smaller than fvy. 



